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Abstract 

We apply the framework of tempered fractional calculus to investigate the 
spatial dispersion of elastic waves in a one-dimensional elastic bar characterized 
by range-dependent nonlocal interactions. The measure of the interaction is 
given by the attenuation kernel present in the constitutive stress-strain relation 
of the bar, which follows from the Kroner-Eringen’s model of nonlocal elasticity. 

We employ a fractional power-law attenuation kernel and spatially temper it, 
to make the model physically valid and mathematically consistent. The spatial 
dispersion relation is derived, but it turns out to be difficult to solve, both 
analytically and numerically. Consequently, we use numerical techniques to 
extract the real and imaginary parts of the complex wavenumber for a wide 
range of frequency values. From the dispersion plots, it is found that the phase 
velocity dispersion of elastic waves in the tempered nonlocal elastic bar is similar 
to that from the time-fractional Zener model. Further, we also examine the 
unusual attenuation pattern obtained for the elastic wave propagation in the 
bar. 
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1 Introduction 

The frequency dependency of the phase velocity and wave attenuation originating 
from the wave-medium interaction is mathematically expressed by the dispersion 
relation k (w), where k and oj are the wavenumber and angular frequency of the wave 
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respectively. Since it is possible to extract the viscoelastic parameters of a material 
from its dispersion curve, the investigation to obtain an accurate dispersion relation 
becomes pivotal for many applications in the field of acoustics, medical imaging, 
seismology and geophysics 

122 02 Sami 122123128]. 

One of the methods which is followed to model dispersion assumes that the ma¬ 
terial possesses memory, which can be classified into two types; memory in time and 
memory in space mnsnM. In materials with time nonlocality or simply mem¬ 
ory, the effect at a given point in space at a given time is dependent on the causes at 
the same point for all preceding times. On the other hand, in materials with space 
nonlocality or spatial memory, the effect at a given point in space at a given time 
is dependent on the causes at all points in space. In this paper, we differentiate 
the dispersion arising due to time nonlocality and spatial nonlocality by referring to 
them as ’’temporal dispersion” and ’’spatial dispersion” respectively. 

For an elastic material with temporal memory, physically it implies that the 
material remembers its past deformations. Such materials do not obey the simple 
Hookean elasticity, but rather they fall into the category of viscoelasticity. The 
classical viscoelastic models, e.g. the Maxwell, Kelvin-Voigt and Zener constitutive 
stress-strain models, lead to temporal dispersion and attenuation of waves in materi¬ 
als mmm- The basic building blocks of the models are elastic springs and viscous 
dashpots, which are arranged in series and (or) parallel combinations to depict the 
interplay between the elastic and viscous properties of a material. The fractional or¬ 
der counterpart of the classical viscoelastic models which are studied using fractional 
calculus, have been found even more useful in investigating the dispersive properties 
of a wide range of complex materials (e.g. biological tissues, polymers and earth 

sediments) 03 [23 123 H3 E3 EH E3 • 

For the last fifty years [23, H) |T?J, it has been well established that under the 
action of a deforming stress, the constituent points of an elastic material display 
range-dependent nonlocal interaction between themselves. Consequently, the applied 
stress is not confined to a local point, but rather distributed to all the interacting 
points of the material. A wave travelling in such a nonlocal material undergoes 
spatial dispersion, in addition to the regular temporal dispersion m • It is essential 
to emphasize that temporal dispersion manifests itself as a result of the localized 
wave-material interaction which is primarily dependent on the frequency to. But 
since space nonlocality in a material is dictated by the wavelength of the traversing 
wave, the resulting spatial dispersion is determined by the wavenumber k. 

Lately, spatial nonlocal operators of the form of fractional Laplacians have been 
used to model dispersive properties of acoustic media mm- However, the study 
of spatial dispersion has not received as much attention as the temporal dispersion 
|16., 18] • The main factor behind this is the satisfactory explanation of dispersive 
properties of most materials by temporal dispersion, which is modelled under the 
framework of classical theory of elasticity. However since the classical theory is 
based on the principles of continuum mechanics, it assumes localized stress-strain in 
the material, and is therefore not applicable in problems where nonlocal interactions 
play a dominant role m- 

The primary motivation for our work is that much of the current research (see, 
e.g. [321 l8l fl8l 13811331 132] 1 has focused more on the mechanical properties of nonlocal 
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elastic materials than its dispersive properties. Second, the results from such studies 
may even influence how new materials can be artificially engineered [3]. Besides, 
we will also illustrate how the framework of tempered fractional calculus can be 
used to overcome the boundary value problems which are often encountered when 
nonlocal elasticity is modelled using fractional calculus. The purpose of this paper is 
therefore to modify the nonlocal elasticity in an infinitely long one-dimensional bar 
with tempering, and then investigate the spatial dispersion of elastic waves in the 
bar. Further, we show how different numerical techniques can be used to solve the 
spatial dispersion relation which is often analytically difficult. 

The rest of the paper is organized as follows. In Section [2l the suitability of 
fractional calculus to study nonlocal elasticity problems is analysed in the light of 
some recent studies. In Section [51 we spatially temper the conventional power-law 
attenuation kernel of the nonlocal elastic bar and examine its significance. Then 
in Section 01 we utilize the mathematical framework of tempered fractional calculus 
to derive the spatial dispersion relation from the constitutive equation of the bar. 
In Section [5J we employ numerical techniques to solve the dispersion equation for 
complex wavenumbers k and subsequently obtain the dispersion plots. Further, we 
compare temporal dispersion and spatial dispersion by pointing out the resemblance 
of the phase velocity dispersion curve obtained from the space-fractional tempered 
nonlocal elastic model with that from the time-fractional Zener model. Moreover, the 
spatial dispersion relation anticipates an unusual attenuation behaviour for the elastic 
wave propagation in tempered nonlocal elastic bar, which we justify by physical 
arguments. Finally in Section [HI the potential implications of this work are used to 
draw some conclusions. 

2 Nonlocal elasticity and the framework of fractional 
calculus 

Nonlocal continuum held theories describe the physics of materials whose behaviour 
at a given local point is determined by the state of all points of the body (see, 
mmmm)- The theory is built upon two assumptions. First, the molecular 
interactions in a material are inherently nonlocal. Second, the theory assumes the 
energy balance law to be valid globally for the entire body. The comprehensive and 
robust formulation of nonlocal theory is ascertained by its wide range of applications 
[16J . Lately, nonlocal theory has also found applications in nanomechanics which 
integrates solid mechanics with atomistic simulations to analyse the strength and 
dispersive properties of carbon nanotubes (CNTs) [33], [30] . 

According to the Kroner-Eringen’s model of nonlocal elasticity, the total strain 
response of a material to a deforming stress a ( x ) is the sum of the local strain e (x) 
and, the nonlocal strain m , m- The nonlocal strain is given by the convolution 
integral of the strain in the neighbourhood e (y) with the attenuation kernel g (x — y), 
where y is the Euclidean distance from the local point x. For an infinitely long, one- 
dinrensional, isotropic, linear nonlocal elastic bar, the constitutive equation is then 
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given as 


cr (x) = E 0 



+ K 


oo 

J e{y)g(x-y)dy , 


Local strain 


Nonlocal strain 


( 2 . 1 ) 


where E$ is the elastic modulus at zero frequency. The material constant k represents 
the strength of nonlocality in the material. As seen, in the limit as k —>• 0, equation 
m reduces to the case of Hooke’s law. It is reasonable to assume the bar to be 
one-dimensional provided the bar thickness on both sides is much smaller than the 
wavelength of the propagating wave. Besides, an infinite length of the bar ensures 
complete attenuation of the wave by the time it reaches its free ends, which makes 
the wave propagation free from reflections and the resulting complexities of bound¬ 
ary value problems. An illustration of the equivalent mechanical model of nonlocal 
elasticity can be found in [521IKl 155], where points of a material are shown connected 
to each other by means of springs of different stiffnesses. 

Although some suitable forms of the attenuation kernel have been suggested in 
the relevant literature (see, e.g. m Ei m suss]), recently there has been a growing 
interest in the power-law kernel (see, [521121 |S1 [551 IT]), which is given as 


a(y) 


l 

r(l-«)|!/r 


( 2 . 2 ) 


where a € [0,1] is the fractional exponent and T (z) is the Euler’s Gamma function 
defined for a complex variable 2 as: 


r 


oo 



x z l e X dx, 9ft (z) > 0. 


(2.3) 


The advantages of opting for a fractional power-law attenuation kernel are three 
fold. First, the exponent a allows a wide range of distance decaying interactions 
to be taken into account and it also plays the role of a scale parameter which is 
inherent to the theory of nonlocal elasticity [52]. Second, the similarity between the 
attenuation kernel and the power-law memory kernel of fractional calculus allows the 
mapping of nonlocal problem under the comprehensive framework of fractional calcu¬ 
lus. Third, power-law solutions are often encountered in experimental measurements 
which further justifies fractional modelling of the nonlocal elasticity. 

Here, it is essential to clarify that the attenuation kernel g (y) should not be con¬ 
fused with the term ” memory kernel” which is often used by the acoustic community 
for the temporal power-law functions. 

Following the commutativity property of convolution and also considering the fact 
that the attenuation kernel is symmetric across the local point of an isotropic bar, 
we substitute equation (EH) in equation m and obtain the constitutive relation 
as, 

OO 


a (x) = E 0 



e(x-y) 

T (1 - a) y a V 


(2.4) 
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Clearly, the absence of a Dirac-delta form of the attenuation kernel followed by the 
presence of an integral in equation (12.411 enforces strong nonlocal elasticity in the 
material. In addition, it is seen that the integral term in equation (12.41) is similar 
to a one-dimensional fractional Laplacian of order (1 — a), which also implies that 
the attenuation kernel given by equation (12.21) has indeed the form of a Riesz kernel 

m H3I- 

Although the framework of fractional calculus offers a powerful tool for inves¬ 
tigating problems characterized by fractional power-law kernels, it suffers with an 
inherent problem of singularity which has also been recently reported by one of its 
proponents in [6j. In the limit as y —>• 0 in equation (12.21) . g (y) —»• oo; this implies 
that the nonlocal contribution to the stress at the local point x is infinite, conse¬ 
quently the displacement at the local point u (0) also grows to infinity, which is 
clearly non-physical. This is evident mathematically as well, since the nonlocal in¬ 
tegral in equation (12.41) cannot pass the p-test of convergence and therefore never 
converges to a real value (see, Chapter 1 in [25]). 

Most works on the topic have encountered the same problem. A recent work 
(see, [32 ] ) which investigates the strength of a nonlocal elastic bar of finite length L 
characterized by the same kernel as given by (12.21) . considers the divergent nature 
of the fractional integral as an open problem. The authors further trace this to the 
original work by Kroner and suggest enforcing the condition: u (0) = u (L) = 0, to 
circumvent the boundary value problem [23|. Further, a very recent article (see, EH) 
overcomes one of the boundary value problems by assuming the nonlocal bar to be of 
infinite length, however it overlooks the divergent behaviour of the nonlocal integral 
at the local point. Another paper on the topic (see, [8]) transforms the problem to the 
case of 1 < a < 2 by expressing the Riesz fractional derivatives in the Marchaud-like 
form. An entirely different solution is suggested in [31], where the author considers 
the distance y to be geodetical, instead of the usual Euclidean. It is worthwhile to 
mention that the article [2] (see, Lemma 1) suggests the possibility of a solution of 
the wave equation for a nonlocal elastic bar in a tempered space distribution which 
does not requires enforcing of boundary conditions. 

To summarise, we find that for the case 0 < a < 1, if homogeneous boundary 
conditions are not imposed, the nonlocal elastic model characterised by the con¬ 
ventional power-law kernel becomes non-physical and mathematically inconsistent. 
Consequently, a closed form expression of the spatial dispersion relation has never 
been obtained for a nonlocal elastic material, and thus the frequency dependency of 
phase velocity and wave attenuation has remained unanalysed. 


3 Spatial tempering of the attenuation kernel 


Now, we follow a similar approach as mentioned in [2] to overcome the boundary 
value problem associated with the fractional nonlocal integral. The method is in¬ 
spired by the exponential truncation of arbitrary large flights produced by Levy 
stable distributions (see, 0126]). Accordingly, we temper the conventional fractional 
attenuation kernel (12.21) as 


1 _ e -%l 

9 (l/) i■ "T"j 7a 

T(l-a)|y| 


(3.1) 
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where a E [0,1] and A > 0 is the tempering parameter. 

As illustrated in Fig. [TJ the tempered attenuation kernel overcomes the sin¬ 
gularity and ensures finite stress-strain at the local point. Besides, the behaviour 
of the tempered kernel in the far nonlocal region is almost identical to that of the 
conventional kernel for all values of a. In the limit as a —>• 0, the applied stress is felt 
all over the bar with little attenuation. On the other hand, in the limit as a —» 1, 
the attenuation kernel approaches the Dirac-delta form implying that the stress is 
completely attenuated in the immediate neighbourhood of the local point. 

As seen, the tempered attenuation kernel does not decay monotonically from 
the local point and therefore, seems contrary to the behaviour expected from an 
ideal attenuation kernel. However we stress that, since the two peaks arising due 
to tempering lie in the immediate vicinity of the local point, the behaviour of the 
material is not affected significantly. Besides, this undesired tempering effect can 
easily be compensated by giving appropriate weighting to the material constant k. 

Replacing the conventional power-law attenuation kernel in equation (12.41) by its 
tempered version m, we then obtain the constitutive relation of the tempered 
nonlocal elastic bar as, 


a (. x ) = E 0 


e(x) + 


2k 


r (1 — a) 


oo 

/ 


(l - e Xy ) e(x - y) 


dy 


(3.2) 


4 Spatial dispersion from tempered fractional calculus 

Since in this section we will transform equation (13.21) to the Fourier domain, it is 
essential to first introduce the mathematical framework of dispersion modelling. For 
a unit amplitude, one-dimensional propagating plane wave, the displacement u is 
given as: 

u(x,t) = e i ( kx ~ ut \ (4.1) 

where x is space and t is time. In a loss-less medium, the phase velocity is c p = ui/k, 
the group velocity is c g = du/dk , and the two are equal. Besides, the wavenumber 
A: is a real quantity, which due to its directional attributes is also called the wave 
propagation vector. As long as the total energy is conserved during the wave-material 
interaction, the convention adopted to model wave attenuation is to assume uj and k 
as real and complex quantities respectively. However, there are also some examples 
where the opposite convention of w as a complex quantity and A: as a real quantity 
have been made (see, e.g. dissipation of thermoelastic waves in solids SI)- In this 
work, we assume the wave propagation vector k to be complex such that 

k = k r + iki, Ay > 0 and ki > 0, (4.2) 

where Ay is the wave velocity vector and ki is the wave attenuation vector. Here, in¬ 
stead of following the common practice of representing the real and imaginary parts 
of k by /3k and atk (see, e.g. [23E3CL3G3EI]), we choose to comply with the nota¬ 
tions mentioned in (14.21) and thereby avoid ambiguities with the fractional exponent 
a of the attenuation kernel. Moreover, the imaginary part ki of the wavenumber A:, 
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Distance, y 


Figure 1: Comparison of the conventional attenuation kernel (top pane) with the 
tempered attenuation kernel (bottom pane) for A = 1. Local point corresponds to 
y = 0. The fractional order a has values 0 (magenta, solid curve), 0.2 (blue, dashed 
curve), 0.5 (black, dotted curve) and 0.8 (red, thicker solid curve). 
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which we denote as ” wave attenuation”, should not be confused with the attenuation 
kernel g(y). Substituting equation (14.21) in equation (14.11) . the displacement of an 
attenuating wave in a dispersive medium becomes 

u(x,t) = e- kiX e i ^ x ~ ut \ (4.3) 


with the phase and group velocities given by the real part of k as c p = ui/k r and 
c g = du/dk r . Both k r and ki are functions of uj and also related to each other by the 
Kramers-Kronig relations due to causality [5Tj. The fact that the Kramers-Kronig 
relations are fundamentally the Hilbert transform pair has two interesting inferences. 
First, both ui and k cannot have imaginary values simultaneously. Second, the wave 
attenuation is always coupled with its velocity. 

Now, we introduce the spatio-temporal Fourier transform as 


CXD oo 

F [u (x,t)] = u (k,uj) = j I u (x,t) e^ kx ^ ut ^dx 


dt. 


(4.4) 


—oo —OO 


As we are concerned with the dependency of wave-material interactions on the wave 
propagation vector k alone, we neglect the time domain and only consider the Fourier 
transform in the spatial domain. Taking the Fourier transform of equation (13.21) . we 


get: 

°(k) = E 0 l^e{k)+ J 2 )|, 

OO 


where 

(4.5) 


h= / F[e(x- y)] y a dy , and 

(4.6) 


0 

oo 



h = j F [e (x — y)} e Xy y a dy. 

(4.7) 


o 

According to the mass conservation principle, strain is given as 


e(x,t) — m = ik ' u ( x ’*)> 

and, its Fourier transform is 

e (k) = ik ■ u (k) . (4.8) 

Next, we consider Newton’s second law of conservation of momentum and the equa¬ 
tion of motion for the plane wave, from which we respectively have: 

d 2 

Vcr (x,t) = P-q^u (x,t) = —oj 2 p u (x,t) , and (4-9) 

d?u (x, t) _ 1 d?u (x, t) 

dF 2 Cg dF ’ 1 j 

where, Co is the phase velocity of the wave at zero frequency and p is the mass density 
of the material. Assuming zero initial conditions and integrating on both sides of 
equation (14.91) . we get 





cr{x,t) = “i 2 jfu{x,t). 

The corresponding Fourier transform of er ( x , t ) is then, 

2 

a (k) = (k). (4.11) 

ik 

Substituting equations (14.111) and (14.81) in (14.51) . we obtain: 

- (k) = E 0 jifc • u (k) + ^ (I\ - I 2 )| ■ (4.12) 

The translational property of the Fourier transform further gives, 

F[e{x-y)\=e~ iky -e(k). (4.13) 

Substituting equations (14.131) and (14.81) in equations (14.61) and (14.71) . we rewrite the 
expressions of I\ and 1 2 as: 


h 


ik ■ u ( k ) 


e- ik yy~ a dy, 


and 


(4-14) 


h 


OO 



e- iky e~ Xy y- a dy. 


(4.15) 


Replacing I± and I 2 in equation (14.51) with their respective equivalent forms (14.141) 
and (14.151) , and then making use of the relation cq = \JEq/ p we readily obtain 



1 + 


2 K 

r(i-a) 



(4.16) 


Next, we proceed to obtain a closed-form solution of the two integrals. It is straight¬ 
forward to see that the two integrals are modified forms of the Gamma function whose 
convergence is established using the dominated convergence theorem (see, Chapter 
3 in [26]). Substituting, x = (A — ik)y in equation (12.31) . where A > 0 and k is the 
complex wavenumber, we obtain a useful expression: 

OO 

n T } Z Lz = [ e~ iky e~ Xy y z ~ 1 dy. (4.17) 

(A + ik) J 
0 

The integrals in equation (14.161) are then evaluated as 

OO 

J e~ iky y- a dy = (ifc) Q_1 r (1 - a ), and (4.18) 

0 
00 

J e- iky e~ Xy y- a dy = (A + ifc) Q_1 T (1 — a ). (4.19) 

0 
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The calculation of similar integrals but for a different exponent of the integrating 
variable can also be found in a recent paper (see, Theorem 2.1 in [36]). Now, substi¬ 
tuting the value of the integrals from equations (14.181) and (14.191) in equation (14.161) , 
we finally have 


U! 

k 2 


— Cn 


l + 2ftj(ifc) Q 1 — (A + ik) a 


(4.20) 


Equation (14.201) is the required spatial dispersion relation which governs the propa¬ 
gation of an elastic wave in a tempered nonlocal elastic bar. It is essential to note 
that the tempering parameter A has a dimension of [L] _1 , one the other hand, the 
material constant ft has a variable dimension of [L] 1_ “. 

Further, in the absence of nonlocality i.e., if A = 0, the spatial dispersion relation 
(14.201) reduces to the classical case of loss-less dispersion given as co = c^k. The 
equivalent form of the constitutive relation a (x) = E^e ( x ) is obtained if the same 
substitution is made in equation (13.21) . 


5 Numerical solution of the dispersion equation and dis¬ 
cussion 

The dispersion equation (14.201) is a nonlinear equation with complex coefficients and 
fractional-order exponents of complex wavenumber k. As it is difficult to obtain the 
closed-form expressions for real and imaginary parts of k, we resort to root finding 
algorithms. Besides, the choice of numerical algorithms is limited since most of the 
available methods are usually applicable for equations containing real coefficients 
and integer-order exponents (see, Chapter 9 in }35j). We even tried to formulate and 
numerically implement the classical root-finding algorithms (such as, the bisection 
method and the Newton-Raphson method) in the complex domain. However, the 
only algorithm which ensured convergence and achieved the desired level of accuracy 
was the Muller’s method. 

Based on the generalization of the secant method, Muller’s method uses quadratic 
interpolation to approximate the given function and then the root of the func¬ 
tion is approximated by the root of the interpolating quadratic (see, Chapter 7 
in m)- The algorithm requires three distinct guesses which are used for three 
functional evaluations to start with, but continues with one function evaluation af¬ 
terwards. The method does not requires evaluation of derivatives of the function 
and its rate of convergence is about 1.84, i.e. nearly quadratic so that the num¬ 
ber of correct decimal places almost doubles with each iteration. Further, we have 
also eliminated the potential round-off errors due to subtractive cancellation by im¬ 
plementing double precision calculation followed by an alternative formulation of 
the roots. For a quadratic equation of the form ax 2 + bx + c = 0, where o, b 
and c are its coefficients, the formula for the root x which we have implemented 
in the algorithm is = —2c/ ^ —\/6 2 — 4ac^, instead of its conventional form 

x\ 2 = b± \/b 2 — 4ac^ / (2a) (see, Chapter 3 in [12]). 

We set the parameters A = 1, Co = 1 and ft = 1 to solve the dispersion equation 
for three different values of a; a = 0.2, 0.5 and 0.8. For each a value, the dispersion 
equation is solved for a wide range of frequency values, 10 -4 < u < 10'. Each 
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frequency decade is uniformly sampled into a minimum of four points, however in 
situations where significant jumps in the values of the roots are observed, a finer 
sampling of six to eight points is followed. The values of the numerically obtained 
roots are accurate up to a minimum of six significant figures. However there was 
an exception for a = 0.8, where the algorithm could not converge properly in the 
frequency range of 10 _1 -10° Hz, which further illustrates the numerical difficulty 
associated with the extraction of complex roots. It should be noted that since the 
coefficients of wavenumber k in equation (14.2011 are complex, roots do not come in 
conjugate pairs. 

As illustrated in Fig. [2] the obtained roots for a given frequency value give the 
respective phase velocity and wave attenuation. The main observations from the 
plots can be summarized into four points. First, the pattern of an increased phase 
velocity with frequency corresponds to anomalous dispersion which is also seen in 
the case of time-fractional Kelvin-Voigt and Zener models (see, Figures 1 and 2 in 
[T9]). Second, in the high frequency regime, the phase velocity levels off just like the 
fractional Zener model. Third, in the low frequency regime of ui < 1CU 1 Hz, wave 
attenuation is relatively higher than in the intermediate-to-high frequency regime. 
Fourth, for frequency values u > 10° Hz, the straight lines in the log-log attenuation 
plot, clearly indicate the power-law dependency of wave attenuation ki on frequency. 
Further, wave attenuation increases as the value of a is increased. 

The obtained dispersion plots can be understood from its underlying physics. 
Unlike the case of temporal dispersion, spatial dispersion is determined by the wave 
propagation vector k of the propagating wave. In the case of low frequency waves, 
the wavelength of the wave is large, as a result the wave interacts with an equivalent 
larger region of the nonlocal elastic bar. The nonlocal attenuation mechanism which 
is spread over the length of the bar leads to greater wave attenuation and slows it 
down. As the frequency of the wave is increased, the size of the wave and hence, 
the nonlocal region of the bar with which the waves interacts becomes smaller. For 
very high frequencies, the wave ’’feels” the material only at a local point and since 
we have only considered the nonlocal attenuation mechanism, wave attenuation ef¬ 
fectively disappears. As explained, a high frequency wave suffers less opposition in 
the nonlocal elastic bar, and therefore traverses with maximum phase velocity and 
little attenuation. 


6 Conclusion 

One of the goals which has been achieved in the present work is the understanding 
of how the phase velocity and wave attenuation are affected as it propagates in a 
nonlocal elastic material. Even though the wave attenuation in a nonlocal bar ap¬ 
pears unusual, it seems very physical. Considering the fact that wave attenuation 
in a nonlocal bar is negligible in the high frequency regime, such a material if en¬ 
gineered could find applications as an effective channel to transfer energy. We have 
also established the importance of tempering and the framework of fractional cal¬ 
culus in investigating nonlocal problems where conventional power-law kernels give 
non-physical results. Further work should include an investigation of the physical im¬ 
plications of the material constant k and the fractional exponent a of the attenuation 
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Figure 2: Frequency-dependent phase velocity (top plane) and wave attenuation 
(bottom plane) in a tempered nonlocal elastic bar. The markers; square for a = 0.2 
(blue, dashed curve), cross for a = 0.5 (black, dotted curve) and circle for a = 0.8 
(red, solid curve), represent the numerically obtained roots. Lack of markers in the 
frequency decade of 10 _1 -10° Hz for a = 0.8 is due to the non-convergence of Muller’s 
algorithm. 
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kernel g (y). 
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